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^vq , Abstract. Three-loop propagator diagrams in HQET can be reduced, using integration by parts, to 

8 basis integrals: 5 trivial ones, two expressible via 3F2, and one only known up to e". Calculation of 
the heavy-quark propagator in HQET is considerably simplified by the non-abelian exponentiation 
theorem. 
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—I. I Heavy Quark Effective Theory (HQET) is an effective field theory approximating QCD 

(N I for problems with a single heavy quark when characteristic momenta are much less 

than its mass m, see the textbook [1]. A method of calculation of two-loop propagator 

diagrams in HQET based on integration by parts [2] has been constructed in [3]. It has 

been extended to three loops in [4]. Methods of perturbative calculations in HQET are 

^ ■ reviewed in [5]. 

^! All generic topologies of three-loop propagator diagrams in HQET are presented in 

Q^l Fig. 1. The corresponding Feynman integrals depend on indices of the lines — the 

degrees of the corresponding denominators. They can also contain powers of numerators 

which cannot be expressed via the denominators. If index of a line becomes zero, 

this line shrinks to a point. In some cases, the resulting diagram contains lower-loop 

/\ • propagator insertions, and is easy to calculate. An algorithm which reduces any three- 

j^ . loop propagator integral in HQET to a linear combination of 8 basis integrals (with 

coefficients being rational functions of space-time dimension d) has been constructed 

and implemented in [4]. The basis integrals are presented in Fig. 2; indices of all lines 

are equal to 1 here. 

The first 5 basis integrals are trivial: 

5i=/3, 52=/l/2, 53=/3, B4=h-^, ^5=/3§, (D 
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are the n-loop sunset HQET and massless integrals. 
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FIGURE 1. Topologies of three-loop propagator diagrams in HQET 



The next two basis integrals are 

56 = Gi/(l,l,l,l,£), 57=/i7(l,l,-l + 2£,l,i; 



(3) 



We use recurrence relations for the two-loop integrals / and J to re-express them via 
convergent integrals (the first two diagrams in Fig. 3) 7(1, 1, 1, 1,2 + e) and 7(1, 1,2 + 
2£, 1, 1). At a non-integer n, the integral 7(1, 1, 1, l,n) has been calculated in [6], and 
7( 1 , 1 , n, 1 , 1 ) — in [4] . In particular, 
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FIGURE 2. Basis integrals 



Using the methods of [7], it is not difficult to obtain the expansions 
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l,2-2£,3 + 4£ 
3 + £,4 + 4£ 



3+(l8C2-y )e + 3(-48C3-7C2 + 75)£^ 



+ 3(216^4 + 56^3 + 1 15^2 - 506)£3 

+ 12(-375^5 + 90^2^3 - 63^4 - 230^3 - 135^2 + V89)£4 + • 



3^2 



l,2-2£,l+4£ 

3 + 2£,2 + 4£ 



= 2 + 6(-2^2 + 3)£+ 12(10^3- 11^2 + 6)£2 

+ 24(-28^4 + 55^3 - 27^2 + 9)e^ 

+ 48(94^5 - 16^2^3 - 154^4 + 135^3 - 45^2 + 12)£^ + • • • 

They can be extended to £^, should it be needed. 



(5) 






2 2 2 
FIGURE 3. Convergent integrals related to Be. ..8 

The most difficult basis integral Bs = 7c (1, 1, 1, 1, 1, 1, 1) can be re-expressed via a 
convergent integral 7c(2,2,2, 1,1,1,1) (the last diagram in Fig. 3). It can be calculated 
in 4-dimensional euclidean coordinate space in cylindrical coordinates. Taking first the 
integral in the radial coordinate of the light vertex, then the integral in its time, and then 
the integrals in the two times of heavy-light vertices, we obtain 



7,(2,2,2, 1,1, 1,1) = 2(^3 -l) + 0(£) 



(6) 



HEAVY-QUARK PROPAGATOR IN HQET 

Due to the non-abelian exponentiation theorem [8], the unrenormalized HQET quark 
propagator in the coordinate space can be written as 



Soit) = -/0(Oexp 

.6 /,v\6e 



'^-(4l^(l)"^^^^(4(l)*<^-^^'-™> 



+ Cf " ^^ f - j {clsAA + CFTpnisif + CATFnisiA + {TFUif snj + • ■ • .(7) 

Suppose we have calculated the one-loop correction to the HQET propagator in the 
coordinate space (Fig. Ad). Let's multiply this correction by itself (Fig. 5). We get an 
integral in t\, t2, t[, t'j with < ti < t2 < t, < t[ < tj < t. Ordering of primed and 
non-primed integration times can be arbitrary. The integration area is subdivided into 
6 regions, corresponding to 6 diagrams in Fig. 5. If the colour factors of the diagrams 
Fig. Ac, d were the same as that of the one-particle-reducible diagram (Fig. Ab), i. e. 
equal to the square of the colour factor Cf of the one-loop diagram (Fig. Aa) (as in the 
abelian case), then the sum of the diagrams Fig. Ab, c, d would be equal to j of the 
square of the one-loop correction (Fig. Aa), as given by the square of the first term in 
the expansion of the exponent (7). In the non-abelian case, however, the colour factor 
of Fig. Ac differs from Cp by —CfCA/2, which is the colour factor of Fig. Ae. This is 
because when we reduce the colour factor of Fig. Ac to that of the reducible diagram 
Fig. Ab, we get an extra term with the commutator [t", t^], which has the colour structure 
of Fig. Ae. Therefore, we should include the contribution of Fig. Ac with —CfCa/1 
instead of its full colour factor into the term sa in the exponent (7). This part of the 
colour factor is called maximally non-abelian or colour-connected [8]. Of course, the 
diagram Fig. Ae also contributes to sa', the diagram Fig. Af, with the one-loop gluon 
self-energy correction, contributes to si (quark loop) and sa (gluon and ghost loops). 
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FIGURE 4. One- and two-loop diagrams for the HQET propagator 

In the same way, we can consider the set of three-loop diagrams obtained by multi- 
plying the corrections of Fig. Aa and /. We can imagine that this set is obtained from 
the one-particle-reducible diagram by allowing the gluon - heavy-quark vertices to slide 
along the heavy-quark line, crossing each other. These diagrams are said to contain two 
connected webs (or c-webs) [8]. Everything is already accounted for by the product of 
the one-loop correction and the (Fig. Af part of) two-loop correction in the expansion 
of the exponent (7), except the contribution of Fig. 6a (and its mirror-symmetric), taken 
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FIGURE 5. Exponentiation theorem 



with the maximally non-abelian part of its colour factor. It contributes to the three-loop 
correction in the exponent. Similarly, out of all the diagrams with two connected webs 
Fig. Aa and e, only those of Fig. 6b, c (plus their mirror-symmetric ones), and d con- 
tribute to saa, with the maximally non-abelian part of their colour factors. This part 
appears, in the case of Fig. 6b, for example, when we commute t" matrices to obtain the 
colour structure of the reducible diagram; it is identical to the colour factor of Fig. 6h, 
equal to CfC\/A. 



^ 








b 




^ 



^--Vv/^ ^-VvT-T^ 



f 




^ 







^ 




J k I m 

FIGURE 6. Three-loop diagrams contributing to the three-loop term in the exponent (7) 



Now we are going to consider all the diagrams with three gluons, both ends of which 
are attached to the heavy-quark line. They are said to contain three c-webs, each of them 
is that of Fig. 4a. We decompose the colour factors of these 15 diagrams in the following 
way. We move the vertices along the heavy-quark line in such a way as to disentangle 
those c-webs. While doing so, we get extra terms from the commutators, having colour 
structures of the corresponding diagrams with the three-gluon vertex. These diagrams 
have fewer c-webs, which are more complicated. Finally, each colour factor can be 
expressed as a linear combination of three ones: Cp (3 c-webs of Fig. Aa); — C^Ca/2 (2 
c-webs of Fig. 4a and e); CpC^/A (1 c-web of Fig. 6h). The first one occurs with the unit 
coefficient in all 15 colour factors. The sum of the corresponding contributions is just 
the term with the cube of the one-loop correction in the expansion of the exponent (7). 
The second colour structure occurs in the diagrams obtained by multiplying Fig. 4a and 
c; the sum of the corresponding contributions is contained in the product of the one-loop 
term and the two-loop one in the expansion of the exponent. We are left with the colour- 
connected parts of the colour factors (a single c-web contributions). They are present in 
the diagrams Fig. 6e, f, g (and its mirror- symmetric). They contribute to saa in (7). 

We are left with the diagrams containing only a single connected web. Those of 
Fig. 6h and i have equal colour factors (this is evident if we close the quark line), 
they contribute to saa- The colour factor of Fig. 6j is zero. This becomes clear if we 
close the quark line and write a three-gluon vertex as the commutator (Fig. 7). The 
diagram with the four-gluon vertex (Fig. 6k) can be decomposed into three terms, with 
colour factors of Fig. 6h, i, j. The diagram Fig. 6/ contains two-loop gluon self-energy 
corrections, including one-particle-reducible ones; it contributes to saa, sif, sia, su. The 
diagram Fig. 6m contains one-loop corrections to the three-gluon vertex, including one- 
particle-reducible ones (i. e., one-loop self-energy corrections to each gluon propagator 
in Fig. 4e); it contributes to saa, sia- 






FIGURE 7. Vanishing colour factor 

In order to find the heavy-quark field renormalization constant, we should re-express 
the unrenormalized propagator (7) in terms of the renormalized coupling as(/j) and 
gauge-fixing parameter a(/j) instead of ^q and ao, and require that the renormalized 
propagator S{t) = Zq So{t) is finite in the limit £ — > 0. After this re-expression, So{t) 
still has the exponential form with the same colour structures as in (7). Therefore, the 
renormalization constant can be written in the exponential form, too: 



Zq = exp 
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The coefficients z are obtained simply by singling out poles in £ from the corresponding 
terms in the exponent for 5o(?) (expressed via the renormalized quantities). 

The heavy-quark field anomalous dimension in HQET is defined as Yg = 
dlogZQ/dlog/j. The /j-dependence comes from OLsd-i) and a{fj) in Zq. Using |3(as) 
and YA(cXi) for their differentiation, we obtain 
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The two-loop term contains no C^ part, and the three-loop one contains no Cp and C^Ca 
parts, as a consequence of the non-abelian exponentiation theorem. The two-loop term 
was obtained in [3]. The three-loop term has been obtained, by a completely different 
method, in [9]; our calculation confirms this result. 

I am grateful to D. J. Broadhurst, K. G. Chetyrkin and T. van Ritbergen for useful 
discussions. 
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